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SUMMARY 


This paper presents preliminary results of a comprehensive study of the structural aspects of a 
numerically simulated number turbulent boundary layer. A direct Navier-Stokes simulation of a flat- 
plate, zero pressure gradient boundary layer at Reg = 670 was used. Most of the known nonrandom, 
coherent features of turbulent boundary layers are confirmed in the simulation, and several new aspects 
of their spatial character are reported. The spatial relationships between many of the various structures 
are described, forming the basis for a more complete kinematical picture of boundary layer physics than 
has been previously known. In particular, the importance of vortex structures of various forms to the 
generation of Reynolds shear stress is investigated. 


INTRODUCTION 


In spite of decades of experimental research, the structure of turbulent boundary layers, and the 
dynamical processes by which turbulence is created and maintained in boundary layers are only partially 
understood. The incomplete nature of our understanding has prevented the results of turbulence structure 
research from making any significant contributions to either predictive modeling schemes or to turbu- 
lence control methodologies. 

This continuing situation has given rise to a novel approach to the problem: a community-wide 
re-evaluation of the knowledge. In simple terms, the aim is to develop a clear picture of what we know 
with certainty and to form a “strategic plan” for attacking the voids in our knowledge. The eventual goal 
is to lay the foundation for a class of statistical-structural” models for turbulent boundary layers. The 
new and critical element of this effort is the involvement of the research community. A detailed 
description of the process and status of the cooperative evaluation project is given in Kline and Robinson 


The limitations of experimental methods for unravelling the puzzle of turbulence physics are well 
understood but are, to a large degree, inherent. Fixed probes at a small number of points in space can be 
weak tools for determining the characteristics of three-dimensional structures that exhibit large 
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variations m time and space. Although many groups agree that simultaneous use of flow-visualization 
and probe arrays are optimum, it remains difficult to resolve Eulerian probe-based quantitative data with 
visual information from markers in a three-dimensional unsteady flow-field. It is to the credit of the 
researchers m the field that so much has been learned about the structure of turbulence through 
expenmentation. s 


To complement the large volume of experimental results in the field, this companion study of the 
cooperative project is designed to take advantage of the direct Navier-Stokes turbulence simulation 
databases available at NASA Ames. Although the simulations are currently limited to simple geometries 
at low Reynolds numbers, there are many advantages to studying turbulence structure numerically as 
discussed later m this paper, and also by Moin (1984). ’ 

Community involvement also plays a fundamental role in the study of the numerical simulation data. 
As the major unknowns are identified through the cooperative evaluation work, their resolution can be 
pursued in the simulations. The questions and suggestions of the participants in the cooperative project 
have helped pose the critical questions for study via the numerical simulation databases. This paper will 
present an overview of the recent results of this approach. F 

The objective of the current effort is to combine the many known structural elements and events with 
new insight gained from the numerical studies, to form a complete dynamical picture of flat-plate 
boundary layer turbulence, at least for low Reynolds numbers. This picture should include information 
on: (1) the spatial character of each class of structure; (2) the spatial relationships between different 
structures; (3) the generation and evolution of structures; (4) the statistical relevance of each class of 

tuSulence^ 5 descnptlons of the dominant sequences of events that are responsible for maintaining 

An additional goal is to investigate the effects of Reynolds number on boundary layer structural fea- 
res, at least oyer the limited range of Reynolds numbers available in the simulations. A final important 
objective is to simulate a variety of experimental techniques in the databases to improve and possibly 
extend the understanding gained from experimental results. 

In the remainder of this paper, a brief description of Spalart’s boundary layer simulation will be 
given as will an outline of the strategy we are employing to study the database. A preliminary taxonomy 
of structures is presented, and examples of each class are illustrated and described. Finally, several 
samples of spatial relationships between structures will be shown. The results confirm a number of pub- 
lished experimental findings and also present new information on the structure of turbulent boundary 


Since the results of this work are voluminous and are still in progress, conclusions concerning time- 
evolutions of structural features will be deferred to a subsequent paper. The results presented herein are 
meant to serve as an introduction to the character of turbulence structure in the simulated boundary 
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ADVANTAGES/DISADVANTAGES OF THE USE OF FULL SIMULATION DATABASES 
FOR THE STUDY OF TURBULENCE STRUCTURE 


The central advantage of the numerical simulations over experimental investigation is that in the 
databases, we can collect a large sample of structural features first, and then decide what statistics are 
pertinent to compute. In the laboratory, the process is usually reversed. The flow can be made visible, 
but quantitative structural information from probes must be somehow extracted by clever statistical 
analysis of the data. This method is necessarily incomplete, and has proven to be seriously misleading 
m some cases. The averaging and filtering inherent in most statistical techniques can discard phase 
information, disguise true characteristics of individual realizations, create false symmetries, and con- 
taminate ensemble averages by the inadvertent inclusion of extraneous flow elements. The current 
availability of numerical turbulence marks the first time in history that a sufficient sample of each 
structure can be inspected a priori and then quantitative statistics gathered for the significant features. 

In addition to this subtle, but crucial point, the full simulation results offer a number of other more 
obvious, but important advantages: 

1. Selected structural features can be pursued both forward and backward in time, and this process 
can be iteratively refined. 

2. All structures of a given type that exist in the computational domain can be made simultaneously 
visible, allowing evaluation of case-to-case variations. 

3. The spatial relationships between the several types of co-existing structures can be studied. 

4. Structural features that are virtually unmeasurable in the lab can be displayed (e.g., low-pressure 
cores of unsteady vortices). 

5. By combining time evolution with the ability to display different types of structures simultane- 
ously, it becomes possible to examine the spatio-temporal relationships among the various elements of 
the dominant structures. 

6. Rapid turnaround time of structural queries and the availability of specific frozen data volumes 
and time-sequences allow iterative reprocessing. 
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The structural aspects of the direct numerical simulation databases have been studied in demh for 

th R T a t ors Kim ; 198 * ^ and Moin ’ 198 * Si 

Spalart, 1988). Recently, expenmentahsts (including the first two authors of this naoerl have bem,n to 
explore the simulations as “numerical wind tunnels" to confirm and extend to 
help design future expenments. Participation by turbulence structure experimentalists in Ae 1987 
Summer Program at the NASA/Stanford Center for Turbulence 
significant findings from investigations of the numerical databases. 

• ,^ ie ] dir f ct Navier-Stokes simulation work is still in its infancy, and has several limitations tw» 
include high costs m time and dollars (both for the simulation itself and for the later analyses) and 

“J 10 " 1 ; low . Reynolds numbers, simple geometries, and periodic boundary conditions Neverthe 
less, the study of the simple flows now available will not only enlighten us 

but™'! also aid m establishing appropriate procedure for application ihe of ,he 


NUMERICAL SIMULATION SPECIFICATIONS 


Navier - Stokes numerical simulation of a flat-plate turbulent boundary layer (Spalart 
i988) has been chosen as the initial simulation database to explore. The code has been run to statistical’ 

as “ fo -^ 

T th R 9 " 670 a ’ due to the hl S h computational cost of the Re fl = 1410 case At Re* - 670 
each time-step is computed on a 384x288x85 grid, comprising 9.4 million nodes At each node for 

Fulton ifl P i » re - SUIe T* a,l . th f ee . com f’ onents of velocity (and thus vorticity) are available Grid 

» viscous lengths m the sheamwise (x) direction, and 4.3 viscous lengths in the sunwise 

Sigridpdl^ 

a " d l Wall : normal tUmensions of 4900, 25M, Sdl l(XJ > ”s- U ' 

will provide sufficient sample sizes for spatial statistics of the major sectoral fe^ures. 

used repeatedly to illusuate several structural features. Many otter sullies t^ftetgtacMedt 
die overall analysts, however, so the ftgu.es presented here should be conside^xampks S not^ 1c 

feature y Me P tetectabk e i™ ; eahratlons . It , is worth notin S that exan >Ptes of each class^f structural 
earure are detectable in nearly every subvolume studied so far. 
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Spalart (1988) has performed a number of statistical checks of the numerical boundaiy layer, which 
generally show good agreement with available experimental results. These include turbulence intensi- 
ties, skewnesses, flatness factors, Reynolds stress and turbulence kinetic energy budgets, and spectra. 
The most significant known discrepancy is a 5% overprediction of the friction coefficient for the highest 
Reynolds number case. In addition, near-wall fourth-order statistics, especially for v', show a depen- 
dence on grid resolution. It is unclear how significant these may be in terms of the structural behavior of 
the turbulence. Results of our analyses so far indicate that the simulation is structurally indis- 
tinguishable from a physical low Reynolds number boundary layer. 


ANALYSIS STRATEGY 


The problem at hand is to understand the physics of turbulence creation and maintenance in bound- 
ary layers. Spalart’s numerical simulation offers access to over 31 gigabytes of flow-field data The 
computational and graphical power available to study the simulation results makes possible the pursuit 
of hundreds of “interesting” questions. Effective progress toward an understanding of turbulence pro- 
duction therefore demands a well thought-out strategy in order to set work priorities. 

We have adopted the following six-phase strategy for the study of the simulation databases (fig. 3): 

I. (a) Identify structures in the numerical boundary layer that are known to exist in laboratory 

boundary layers. These can be grouped into eight classes, which are delineated in the next 
section. 

(b) Once examples of these structures have been identified in the simulation results, obtain detailed 
descriptions of their spatial character by studying individual time-steps. 

n. Determine the spatial relationships between the various structures by visualizing them together in 
individual “frozen” time steps. 

HI. Study the spatio-temporal relationships between structures by using all the stored time-steps and, 
where necessary, computing intermediate solutions at very small time increments (At + = 0.3) 
between the stored time-steps. 

IV. E>etermine the statistical significance of the various structures, their apparent generation mecha- 
nisms and any interaction events. Populations, frequencies in space and time, and statistics of 
occurrence will be computed. Contributions to future statistical/structural models will be made at 
this stage by computing probability distributions of events and structures, and their contributions to 
specific terms in the Reynolds-averaged model equations. 

V. Combine the results of stages I through V to develop a phenomenological model of the maintenance 
of turbulence in flat-plate boundary layers. This should include a description of the important 
sequences of events, and should describe the roles of all classes of experimentally observed struc- 
ture. The flow of cause and effect within these sequences should be modeled, as well as any sig- 
nificant interactions between structural features. 
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VI. Simulate various experimental techniques, including multi-sensor probes and convected visual 
markers. 

Phases I through IV involve only the kinematics and statistics of turbulence structure. As J. M. 
Wallace has reminded us, a picture of the dynamics does not emerge from passive observation of phe- 
nomena. Some kind of model must be proposed which links observed sequential events with causative 
mechanisms. Thus, only Phase V directly addresses the problem of dynamics in the turbulent boundary 
layer. J 

In the context of this strategy, most quantitative experimental results have addressed Phases I, II, and 
IV, while most flow-visualization studies could be categorized as Phase in. A number of authors have 
proposed structural models which include dynamical hypotheses and so belong in Phase V. Our Phase 
V work will involve testing the many models published by the community in the simulated turbulence, 
with the goal of combining the best parts of each into a single, and hopefully more complete, phe- 
nomenological model. 

Phase VI of the strategy involves the simulation of various probe-based and visual experimental 
methods. In the numerical simulation, spatial variations can be considered separately from temporal 
variations. This may lead to clarification of probe-based results, which are generally Eulerian and hence 
treat space and time as a single, inherently combined variable. 

The six-phase analysis strategy allows us to build up knowledge in a logical order, and facilitates 
prioritization of the large number of possible queries that might otherwise become disorderly and over- 
whelming. In the current study, the phases of the strategy represent an iterative, rather than a one-pass, 
sequential process. So far, work has been done under all six phases, but none have been completed. As 
noted, the current paper focuses largely upon results of Phases I and II of the research program. 


TAXONOMY OF STRUCTURES 


A consistent and accepted definition of “structure” (in the turbulence sense) has not yet surfaced. 
Thus, it may be viewed as a hopeless task to categorize the many flavors of turbulence structure to be 
found in the literature. However, a taxonomy, even though somewhat arbitrary, can help focus thinking 
and inject organization into the somewhat tangled mass of information on the subject. 


Therefore, we have divided the various structures, events, flow modules, and organized motions 
found in the literature into eight classes. The intent is to include all types of non-random events that 
have been singled out for study by the researchers in the field, rather than utilize a definition that would 
necessarily exclude some of the data. This classification is only one of several possibilities, and the fol- 
lowing list is not a conclusion of the current work, but rather a starting point for it. 

A Classifica tion of Boundary Laver Turbulence Structures : 

( 1 ) Wall low- speed streaks 

(2) Ejections of low-speed fluid outward from the wall 
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(3) Sweeps of high-speed fluid inward toward the wall 

(4) Vortical structures of various forms 

(5) Near-wall shear layers, exhibiting strong spanwise vorticity and du'/dx 

(6) Near-wall “pockets,” seen as regions swept clean of near-wall marked fluid in experiments 

(7) Large (delta-scale) discontinuities in the streamwise velocity, or “backs” 

(8) Large (delta-scale) motions capped by bulges in the outer turbulent/potential interface 

Kline and Robinson (1988) give descriptions and some of the known characteristics for each of these 
structural classes. 

The present classification is meant to provide a skeleton of organization on which the study of the 
n u m e n u cal simulations can build. Reflected in the broad coverage of the list is the objective of including 
all of the known structural features of turbulence in the current study, rather than just a few. Thus the 
categorization is comprehensive if not final. 


PHASES I AND H: EXAMPLES OF THE SPATIAL CHARACTER OF EACH STRUCTURAL 
CLASS AND SELECTED SPATIAL RELATIONSHIPS 


Sublayer Streaky Structure 

Figure 5 shows an instantaneous plan view of the x-z plane at y+ = 2 in the simulated turbulent 
boundary layer. The low-speed regions (streamwise velocity less than the local mean) are elongated and 
thin, m agreement with many experimental observations. (Here and elsewhere in the paper, “local 
mean refers to the spatial average over the entire computational x-z plane at constant y + l for at least 
five time steps.) Streamwise coherence of single low-speed streaks often exceeds 1500 viscous lengths 
with widths rangmg from 20 to 80 Az+. The streaky character of the instantaneous low-speed regions is 
somewhat more pronounced than shown by marked fluid in laboratory flows. This is because markers in 
low-speed streaks often lift away from the wall during ejections, leaving the remaining near- wall streak 
unmarked and invisible However, the current results appear consistent with the wall-temperature data 

of nirata and Kasagi (1979) m which continuous visualization at the wall was achieved using liauid 
crystals. b H 

. high-speed regions in the sublayer (fig. 5) are considerably less elongated and somewhat wider 
an the low-speed regions. Streamwise lengths of high-speed regions seldom exceed 800 x + and 
widths range from 40 to 1 10 Az+. This difference in “streakiness” is to be expected, since high-speed 
fluid generally originates from outside the sublayer where the mean velocity gradient is lower and the 
flow is known to be less streak-like. Experimental flow visualizations have not shown the character of 
high-speed regions clearly, however, since fluid markers introduced near the wall collect in the low- 
spee regions. Since high-speed regions in the sublayer are less elongated than the low-speed regions 
the term high-speed streaks” is not used in the current study. 
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Another -characteristic of low-speed streaks is shown in figure 6, which is an instantaneous plan-view 
of the boundary layer at y+ = 15. The colored contours mark regions of significant value of the prod- 
ucts (— u X+w ) and (-u')(-w’). This product serves as a “streak-turning” parameter, since it highlights 
low-speed streaks (-u') that simultaneously exhibit spanwise motion (±w’>. This correlation may be 
termed u ’w'2 and u’w’3, after the quadrants it occupies in the u'w' hodograph plane. It was noticed by 
one of the authors and J. Swearingen that u'w'2 and u'w'3 contours display a vaguely “cross-hatched” 
pattern, suggestive of possible transverse propagation mechanisms near the wall. Swearingen et al. 

£. . . ) have drawn parallels between this behavior and spanwise oscillations related to breakdown insta- 
bilities in curved-wall flows. In the simulated boundary layer, significant near-wall spanwise velocities 
are commonly associated with “leg” or “neck” vortices (see below). A complete explanation for the 
spanwise motions of the low-speed streaks awaits further investigation. 


Ejections 


The quadrant technique will be used to define and identify ejections, or regions of low streamwise 
velocity with motion outward from the wall. The quadrant method involves splitting the instantaneous 
product of u' and v' into four categories according to the quadrant occupied in the u’v' hodograph 
Since its introduction by Wallace et al. (1972) and by Willmarth and Lu (1972), the quadrant method has 
been widely used to assign structural significance to motions that contribute to the boundary layer clo- 
sure term, the Reynolds shear stress. 

It is important to distinguish between the Reynolds shear stress -uV", which is defined only as an 
average, and the instantaneous value of u V. To facilitate this distinction, the new term “shear product” 
will be used to refer to the instantaneous value of u V at a location in space. 

Figure 4 shows the four quadrants into which the shear product can be split, and the names associ- 
ated with their respective motions (Wallace et al., 1972). The contributors to the positive Reynolds 
shear stress are quadrants 2 and 4. The concepts of ejection and sweep imply a coherent movement of 
fluid, rather than the occurrence of quadrant 2 or 4 motion at a single point in space. In the current 
work, ejections are arbitrarily defined as instantaneous regions of quadrant 2 shear product that extend at 
least ten viscous lengths in any direction. Sweeps are defined similarly. Thus, every ejection is a quad- 
rant 2 motion, but not necessarily vice versa. Note that under this definition, sweeps and ejections mav 
occur any where in the boundary layer. J 

Figure 7 shows contour surfaces of quadrant 2 shear product in the instantaneous computational sub- 
volume, S. For the figure, the contour value was set to detect strong quadrant 2 motions, with u’v'2 
equal to 4u^ . This and many other realizations show that strong ejections occur with two types of 
instantaneous spatial character: 1) thin and elongated in the buffer and lower log regions, and 2) broader 
and more blob-like in the outer log and wake regions. The spanwise extent of the outer-region ejections 
ranges up to Az - 200. Since ejections near the wall originate from low-speed streaks by definition 
then- character is expected to be at least somewhat elongated. Ejections are not observed to extend over 

the entire streamwise length of a streak, however. Only portions of the low-speed regions extend out of 
the sublayer. 


Sweeps 


Sweeps are regions in which the shear product uV occupies the fourth quadrant in the u V plane 
(fig. 4). Instantaneous contour surfaces of quadrant 4 shear product in the subvolume S are shown in 
figure 8. The contour threshold is the same as for figure 7; that is, the iso-surfaces correspond to 
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quadrant 4 shear product equal to 4u t 2 . At any given time, strong sweeps occupy less volume in the 
simulated boundary layer than equally strong ejections. This spatial result agrees with the probe-based 
statistics of Wallace i et al. (1972), Willmarth and Lu (1972), and others, who have shown ejections to 
donunate the contributions to the Reynolds shear stress outside the buffer region. The spatial character 
ol high-speed fluid with a wallward velocity component is somewhat different than that of ejections, 
strong sweeps are only occasionally elongated in the streamwise direction, and tend to have span wise 
dimensions of less than 100 Az + , regardless of distance from the wall. 


Low-Speed Streaks with Ejections and Sweeps 

Figure 9 is a plan-view of an x-z plane in the simulated boundary layer at y + = 15. Ejections and 
sweeps are shown for the same contour levels as listed above. Low-speed streaks are identified by con- 
ours of streamwise velocity 3u x below the local mean. Strong ejections are seen to occur only on low- 
speed streaks, as they must by definition (both are regions of -u’). The simulated boundary layer is seen 
to agree with Conno and Brodkey’s (1969) and Bogard and Tiederman’s (1986) finding that several 
ejections may arise from a single low-speed streak, although the number of ejections per streak is 
dependent upon the contour levels chosen for -u’ and uV 2 . The contour value for second quadrant u'v’ 
in figure 9 corresponds to a value of approximately 2.6, when normalized by the product of uYrms) and 

v (rms) at y - 15. Strong ejections seem to originate from the central portion of streaks more often 
than trom either end. 

As expected from experimental results and from previous simulations (Moin, 1984), the contours of 
strong quadrant 2 and 4 motions are highly intermittent in space; neither displays the extended stream- 
wise coherence of the low-speed regions, although quadrant 2 motions are somewhat elongated in the 
x-direction. The spanwise scale of sweeps is slightly larger than the low-speed streaks and the ejections 

than 80 z+ - Th “- swee > in ™ a S2S3 

o,h^wf an "l fl °r' and iS e ^ dem in fisure 9 - Wps and <*«*»• virtually ueTfotow each 
l ^ m the streamwise direction, contrary to several published structural models. However 

lowed rwT g ° f t the i Wee P /eje u Ctl0n xf ir in the x ' z plane wouId cause a fixed probe to detect one fol- 
InL d 1 ^ m 6 by thC °? er - Most sin § le 'P oint conditional sampling schemes do not allow for 

spanwise variations, so ensemble averages that show sweeps following ejections are to be expected but 
are misleading if mteipreted via Taylor’s hypothesis, with the additional assumption that the motion 
lacks significant spanwise component. 

The lateral pairing of vertical motions near the wall is more clearly shown in figure 10 in which 

isrsr v ' x - z pianeat y+ = 15 R * ure io sk-ws 

Moin f logi in^nn^ifl^ ° ccur 1 cl0sely beside inward vertical motions (-v’). This was also noted by 
om (I 984 ) in channel flow simulations. This character suggests that in the buffer region the occur- 
rence of strong outward motion is strongly associated with an occurrence of wallward motion The 
l?nath 1SC ^ imensions of agnrficanl wallward motions are limited to less than approximately 70 viscous 
whkSh a Th K Patt ? n m fl ?T re 10 1S hlghly suggestive of Realized, quasi-streamwise vortical motions 
197 ] h Ri dlSCU , S ?1 1T ] a number o f Publications (Bakewell and Lumley, 1967; H. T. Kim et al’ 

^ na frs mnT, f er 19 ? ; S T ith ** 1983) ‘ The spanwise dimensions of & 

p s range from 40 to 70 viscous lengths. The streamwise extent of the pairs ranges from 50 to 

250 viscous lengths; however, this is not indicative of the streamwise dimensions of near- wall 


9 



quasi-streamwise vortices. Since the figure depicts a plane of data at constant y, a vortex with an 
inclination angle to the wall would be intersected in a limited region by the plane. Note that laterally 
arranged triplets of alternating positive and negative v\ which suggest the occurrence of vortex pairs 
are relatively uncommon, but do exist. Most significant vertical motion in the buffer rcgZ^pZl lo 
be associated with single vortices. 6 appear iu 


Vortical Structures 

For the current work, vortical structures are defined as identifiable vortices, usually with segments of 
varying orientation m space. Under this definition, regions of significant vorticity do not necessarily 
qualify as vortical structures. The problem that arises from any reference to vortices is the definition of 

fav^w UnSteady V1 ?S US flOW ‘ Ir ? the cross ' stream (y- z ) P lane of a zero pressure-gradient boundary 
layer the time-averaged flow is nearly zero. In this plane, a working definition of a vortex is as follows- 
A vortex exists when instantaneous streamlines mapped onto the cross-stream plane display a roughly 
circular or spiral pattern; the core of this vortex will be oriented substantially streamwise. In any plane 
with significant mean motion, this definition becomes ambiguous due to the necessary choice of a refer- 

h m M Ve T lty 'u Fr ° m areference frame movin § with the core of the vortex, this same definition 
should hold: although it requires identifying the vortex core before defining the vortex itself Thus the 
above definition of a vortex is useful as a check of other vortex-identification methods, but not as an a 
pnon tool m any plane other than the cross-stream (y-z) plane. 

Distinction must be made, of course, between vortices and vorticity. Since regions of strong vortic- 
ity are not uniquely associated with vortices, the spatial character of the vorticity field in the turbulent 
boundary layer tends to be noisy and inconclusive with respect to vortical structures. This is demon- 
strated in figure 1 1, m which iso-surfaces of the three fluctuating components of vorticity in subvolume 
. are shown. Clearly, understanding of vortical structures requires a technique for marking only 
instantaneous vortices in the turbulent flow field. 6 y 

Vortex lines have been used extensively by Kim and Moin (1986) to study the structure of vortices 
m a numerically simulated turbulent channel flow. This method, however, requires some meaningful 

W ;i W , hlCh , t0 Ch °°f S I ming locations for the vortex line tracing. In addition, vortex lines do 
not reflect the local strength of vorticity, so that dense bundles that occur in regions of very low vorticity 
could attract undue attention. Finally, although vortex lines must move with the fluid particles (neglect- 
ing viscosity) during time evolution, vortices, which are bunched sections of vortex lines, need not 
thus the method of vortex lines can be arbitrary and inconclusive unless coupled with an independent 
means of identifying vortex locations. 

Near-wall vortices with mostly streamwise orientation have been the subject of a number of experi- 
mental investigations (Praturi and Brodkey, 1978; Blackwelder and Eckelman, 1979; Smith and 

i i Swean "gf and Blackwelder, 1987). By taking a cross-stream cut through the simu- 
ated boundary layer and plotting velocity vectors or streamlines, such vortical motions are easily found 

and iKn^h^ defm iI-T of a vortex - Fl § ure 12 show s an instantaneous cut through a low-speed streak 

t . g . ? nn f hl § h 'Speed region m subvolume S. Velocity vectors in the cross-plane (vV) clearly 
meet the criteria of roughly circular streamlines and so a trial vortex has been identified in the data. 

. na _t ls req J lired to t i race the instantaneous, three-dimensional extent of vortical structures in 

sptice Our proposal is to utilize the fact that vortex cores should be regions of low pressure. In fig- 

dll 13 ’ ? W ~ prCS T e and cross-plane velocity vectors are plotted for the same instantaneous plane of 
data as shown in figure 12. Fluctuating pressure is shown in color in figure 13, with lighter shades 
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indicating lower pressure. The instantaneous v'w' velocity vectors are overlaid as in figure 12. 
Marking the low pressure clearly identifies the vortex core in the example. In figure 13, the radial 
pressure gradient (normal to the core) is fairly large at the outer edge of the core, so the value of the 
pressure contour used to locate the vortex is not critical. A white contour line has been drawn in figure 
13 for values of pressure equal to 4pu x 2 below the local mean. (Local mean refers to the average value 
at each y-location over the entire computational x-z plane, for several time-steps of the simulation.) 
For reference, the rms of the pressure fluctuations, normalized by pu x 2, varies from about 2.5 near the 
wall to 0.5 at the mean boundary layer edge (Spalart, 1988). 

In three dimensions, contour lines become contour surfaces. The 3-D contour surface of constant 
low pressure that corresponds to the 2-D contour line in figure 13 is shown for subvolume S in figure 
14. As can be seen m figure 14, low-pressure regions tend to be elongated and occur in a variety of 
shapes and sizes, apparently indicating a wide variation in the three-dimensional form of vortical struc- 
tures. This observation is confirmed by the study of many other subvolumes in the database. Compared 
to the vorticity field shown in figure 1 1, the instantaneous pressure field of figure 14 is relatively noise- 
free, with a variety of structural features clearly evident. (Figures 1 1 and 14 display the same data set ) 

Comparison of figures 1 1 and 14 emphasizes the importance of distinguishing vortical structures from 
regions of concentrated vorticity. 


The technique of marking low-pressure regions allows application of the vortex definition advanced 
earlier. The velocity of a candidate low-pressure region is first determined. Then the instantaneous 
streamlines in a plane normal to the axis of elongation of the low-pressure zone are computed from a 

frame of reference moving with the low-pressure region. Roughly circular or spiral streamlines confirm 
the existence of a vortex. 

It is not immediately clear whether a discernible low-pressure region is both a necessary and a suffi- 
cient condition for the existence of a vortex in the flow field. Many cuts normal to the axes of elongated 
iow-pressure regions have been made. In every case so far, the velocity vectors in a cut-plane moving 
with the mean velocity of the low-pressure region has displayed roughly circular or spiral instantaneous 
streamlines. Although most low-pressure regions are elongated, there are several convoluted volumes of 
low pressure that are globule-shaped. So far, the analyses show that virtually all elongated low-pressure 
regions are vortices, pie non-elongated regions are not always identifiable as vortices, and are assumed 
to be due to other fluid interactions, such as local accelerations. Thus, the technique of mapping elon- 
gated low-pressure regions provides a means of identifying the instantaneous spatial extent of vortical 
structures in the simulation. 

The loop-like vortical structures shown in figure 14 are strongly suggestive of the hairpin or horse- 
shoe vortices discussed in many previous papers, beginning with Theodorsen (1952). To facilitate dis- 
cussion of the vortical structure topology, various sections of the structures are labeled as follows: heads 
are vortices with mostly transverse orientation; necks are vortices that generally extend wallward and 
upstream from heads at an angle of approximately 45 degrees to the wall; legs are vortices with mostly 
stream wise orientation. Legs sometimes appear connected to and trailing from a head/neck vortical 


though statistical compilations have not yet been completed, some preliminary statements may be 
made concerning the structural nature of the elongated low-pressure regions. First, arch- or horseshoe- 
ke vortical structures (head/neck) combinations are common, but long leg vortices are not. The legs 
that do exist can be quite long in the streamwise direction, reaching up to at least Ax + = 400. True 
airpm vortices, with two trailing legs, are very rare. Many vortical structures are strongly asymmetric. 
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Hook-like shapes, as in figure 14, appear relatively often. In some cases, the missing neck can be 
identified on an asymmetric hook-like vortical structure when the contour level of the isobaric surfaces 
is varied. Legs can be found without necks, and are seen both near the wall (center at y+ < 60) and in 
the wake region. Vortical arches occur with a variety of scales, and their heads are found at y-locations 
ranging from the buffer region to the edge of the boundaiy layer. Arches are generally as tall (in the 
wall-normal direction) as they are wide. The spanwise extent of vortical arches ranges from 50 to 
350 viscous lengths in the simulated layer. The above conclusions hold true over a range of — p 1 contour 
values. In general, varying the -p' contour level affects the diameter and, in some cases, the topologi- 
cal connectivity of the low-pressure regions, but does not significantly alter the overall form of the 
vortical structures marked by elongated regions of low pressure. A variety of statistics on the spatial 
character of vortical structures in the boundary layer will be reported in subsequent publications. 


Vortices with Low-Speed Streaks, Sweeps, and Ejections 

An example of the spatial relationship between low-pressure regions, low-speed regions, and quad- 
rant 2 shear product is shown for the sample subvolume S in figure 15. Low-speed regions are identi- 
fied by iso-surfaces of streamwise velocity 3u^ below the local mean. Surface contour values for uv 2 
and pressure are the same as for figures 7 and 14, respectively. Portions of at least two elongated low- 
speed streaks are visible near the wall, as are significant regions of locally low-speed fluid away from 
the wall, both within and beyond the buffer zone. Figure 15 shows low- speed fluid concentrated below 
and upstream of the vortical heads, and alongside the trailing vortex leg. The arch-like vortical struc- 
tures in the figure straddle elongated low-speed regions, similar to the manner described by Runstadler 
et al. (1963) and by Smith (1984). Just upstream of the smaller of the two vortical arches is a small 
transverse region of low pressure, which appears to be a newly formed secondary vortex, possibly rolled 
up on the low-speed fluid lifted by the action the downstream vortical arch in the manner described by 
Acarlar and Smith (1987). The formation of this particular secondary vortex is clearly visible in time 
evolutions of the simulation. 


The regions where low-speed fluid possesses a significant outward velocity component are marked 
m red in figure 15 as quadrant 2 ejections. Most regions of lifting low-speed fluid are in close proximity 
to the isobaric surface. 

The spatial relationship between low-pressure regions, strong ejections, and strong sweeps is shown 
forsubvolume S in figures 16 (oblique view) and 17 (plan view). The instantaneous spatial association 
of both sweeps and ejections with the vortical structures is striking. 

Strong ejection regions (uV '2 > 4u T 2 in fig. 16) apparently occur in two locations: (1) Alongside 
quasi-streamwise leg vortices, on the inboard (towards the head) side of the leg. This is the upward- 
rotating side of the vortex, so the picture is kinematically consistent. (2) Underneath and upstream of a 
head vortex. Ejections in location 2 are significantly farther from the wall than those in location 1. Both 
the head and neck seem to contribute to the outward induction of low-speed fluid from upstream and 
below the head. Since leg vortices are relatively rarer than head/neck combinations, the majority of the 
strong ejections seem to be associated with the head/neck vortical structure. 

For ejections along quasi-streamwise vortices, figure 17 clearly shows that a single vortex can be 
responsible for the lifting of low-speed fluid, in agreement with Pearson and Abernathy’s (1984) 
numencal studies of vortices in viscous shear layers, and with the ensemble-averaged simulation results 
of Guezennec et al. (1987). This view differs from hairpin vortex models (e.g., Offen and Kline, 1975; 
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Smith, 1984) of the boundary layer, which have fluid lifting due to the combined action of a pair of 
counter-rotating quasi-streamwise vortices. 

Strong sweeps are observed mainly in two locations: (1) On the outboard, wallward-rotating side of 
the leg vortex. Here, high-speed outer fluid is induced toward the wall in a sweep motion. (2) A wall- 
ward induction of high-speed fluid by the outboard side (away from the head) of a neck vortex. Sweeps 
are more common in location 2, since neck vortices are more common than leg vortices in the simulated 

boundary layer. In either case, sweeps are found to be closely related to the existence of vortical 
structures. 


, . £ § e Percentage of the Reynolds shear stress in the boundary layer is known to be produced by 
highly localized strong ejections and sweeps (e.g., H T. Kim et al., 1971; Brodkey et al., 1974- 
Willmarth, 1975). Visualizations of the simulated boundaiy layer (fig. 16 is an example) suggest that 
the generation of Reynolds shear stress is closely tied to the existence of coherent vortical structures 
Ihe conclusion that vortical structures are a critical element of turbulence production and momentum 
transfer in the boundary layer is in agreement with a number of previous experimental and numerical 


Near-Wall Shear Layers 

The existence of sloping du'/dy shear layers near the wall (y+ < 80) of turbulent channels and 
oundary layers has been established both by flow- visualization studies (Corino and Brodkey 1969) and 
by probe-based results (Kreplin and Eckelmann, 1979; Johansson et al. 1987a). The VITA technique in 
particular has been employed in a number of experimental studies to detect rapid changes in the stream- 
wise velocity Correlations of the experimental data have shown these velocity jumps to correspond to 
the passage of local shear layers with a steep downstream slope. Bogard and Tiederman (1987) have 
shown that such shear layers occur on the upstream face of an ejection from the buffer region at the 
interface between the lifting low-speed fluid and the higher-speed fluid overtaking it from behind. The 
limited ability of the VITA technique to detect ejections implies that near-wall shear layers also occur in 

stood however ^ dynamical si S niflcance of these shear layers has remained poorly under- 

Recent work by Jimenez et al. (1987) in numerically simulated channel flows has shed new light on 
the role of near-wall shear layers. These authors have proposed a mechanism for the generation and 
maintenance of the shear layers which is essentially equivalent to that responsible for the instability of 
two-dimensional T ollmien- S chlich ting waves. 

Additionalrecent work on near-wall shear layers in numerically simulated turbulent flows has been 
reported by Johansson et al. (1987b). These authors found that shear layers in the near-wall region 
?rwf §ate Wlt ! 1 a V f l0cl T ty of about 106 U T’ and retain their coherence over streamwise distances of 
™ :- OUS le ?. g . ths - J° hansson et al- also show that near-wall shear layers make important contribu- 
tions to the conditionally averaged production of turbulent kinetic energy. 

In the current study of the simulated boundary layer, the sublayer and buffer region are found to be 
near ' J wa11 shear Shear layers can be found in any randomly selected x y 
slice through the data, and are most common below y+ = 80. Several examples of near-wall shear lay- 

7l\e Ho thfw^nnTh ^ i?° Wn m l hT r Unrdated side - views in 18. The structures are nearly par- 
allel to the wall in the sublayer and slope at an increasing but shallow (<20°) angle farther from the wall 
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The spanwisc dimension of most near-wall shear layers ranges from 30 to 60 viscous lengths imnlvW 
that highly three-dimensional mechanisms are involved in their formation. ’ P y S 

Near- Wall Shear Layers with Vortices and Ejections 

A new observation made in the current study is that some of the near-wall shear layers roll un into 
^ansverse vomces Transverse roll-^ of internal shear layers has previously to 

taS?.™ P T'?“ tS b , y o N l Ch l S “ f • (1973) : Fi « ures 19 md 20 show an example of 
•. A , . i . ^ . . ?■ ^ figure 19, the shear layer is visualized by contours of span wise vortic 

lty. Although the vorticity is high everywhere along the shear layer, a low-pressure region is seen to 

fhethearV 1 thC ° Uter r tip ’ ™ bere a transverse vortex is indicated. The existence of a vortex at the tip of 
the shear layer is confirmed by the instantaneous vector field in figure 20, which has been plotted in a 

reference frame moving downstream with the shear layer tip at its propagation velocity of 15 u T Obser 
vation of computed time-sequences confirms that the particular transverse vortex in figure 20 does 
appear after the formation of the shear layer. ® 

• °/ ^nndrant 2 shear product show ejection-type contributions to the Reynolds shear stress 

UD of die shear , reg101 ? ' . stron g quadrant 2 motion is observed to be associated with the transverse roll- 
h u h u 1 yer ’ but 1S not present alon £ the remainder of the shear layer. Although the dynamic 

fication 8 Ae ^ ^ ^ r ° U ' U ? the ejection ° f low - s P eed streak fluid needs clari- 

1 087 - c’ • U3Tent f ^ 1 1 ng f su PP ort the suggestion of several experimental studies (Acarlar and Smith 

wil se ", and BladtwddCT - Johansson et al.. 1987a) that shear-Iay4 SSte 

wall region may play a major role in the production of turbulence. inxianiuttes in the 

High-Pressure “Potatoes” with Low-Pressure Vortical Structures 

«ti,H?l?f SimU ^ 0n data , base Presents opportunities not available in the laboratory, and this allows the 
study of an additional class of “structure,” namely high-pressure regions. Figure 21 shows die Z ig£- 

subvoIume s of the “ * ada > » 

♦ * nd Iow 'P re ssure regions in the simulated boundary layer exhibit very different characters 
noted by Mom and Spalart (1988). Whereas the low-pressum regions tend to^ 

P ure regions are roundish and potato-shaped. These “pressure potatoes” often occupy the space iust 
upstream of the head portions of arch-like vortical structures, and appear to exist at a converted local 
stagnatmn point which lies at the interface of high- and low-speed fluid. The dynlic ^au^Lf 

h^h nrp«° nvecting P ressure sr adl&nt ret up between the low-pressure vortex core and the following 
high pressure region have yet to be fully explored. ® 

Wall-Pressure Fluctuations 

freure°22° U TT,? “ taneous -Pressure fluctuations are shown for the simulated boundary layer in 

(Wn^efack^et^d^tOTT^and 1111 ^ ” th , an . elon « ated - in agreement with experimental results 
uils a, et al., 1977) and previous stmulanons (Giotzbach and Schumann, 1979: Moin 1984) 

High- and low-pressure regions appear similar in dimensions, with most regions ranging from 50 to 
200 wail units in both directions. However, occasional high- and low-pressure regions Ihh extents of 
over 400 vtscous lengths occur. Although wall-pressure disturbances of this reafe would Ire expS to 
be associated with events in the outer flow, their specific sources are unclear P 
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Pockets 


When distributed markers are introduced into the sublayer of a turbulent boundary layer or channel 
roughly circular regions devoid of marked fluid appear in the plan view. These have been named pock- 
ets (Falco, 1980a), and they give the visual impression of being a “footprint” of some outer structure that 
induces fluid toward the wall. Kim et al. (1987) have shown that pockets become visible in numerically 
simulated channel flow only with distributed fluid markers, not with simulated bubble-lines. This 
observation confirms a similar experimental finding of Falco. However, time evolution of the marked 
flow field is necessary to detect pockets with distributed markers. In the current analysis of the simu- 
lated boundary layer, we have endeavored to detect pockets in instantaneous flow fields. 


If pockets are associated with wallward-moving fluid in the sublayer, the instantaneous streamlines 
in the x-z plane would be expected to diverge in the pocket region. Figure 23 shows instantaneous 
streamlines computed in the y + = 2 plane of the simulated boundary layer. Regions of abrupt spanwise 
divergence of the streamlines are evident in the figure. The spanwise dimensions of these regions are 
y b etween and 100 viscous lengths, in agreement with pocket dimensions reported by Falco 
(1983). A time evolution of simulated fluid markers is being performed by the authors to check whether 
the regions of streamline divergence in the numerical boundary layer are indeed pockets. 


Pockets with Sweeps 

Spanwise divergence of instantaneous streamlines in the sublayer is, by continuity and the streaky 
nature of the streamwise velocity in the sublayer, most likely to be associated with wallward velocity 
Wallward fluid in the sublayer arrives from regions with relatively higher streamwise velocity, so 
diverging instantaneous streamlines are generally expected to be associated with UV 4 motions, or 
sweeps. Since sweeps are observed to be commonly associated with vortices in the simulated boundary 
layer, pockets may be the footprint of fluid swept toward the wall by vortical structures. Preliminary 
studies in the simulated data suggest that this reasoning is correct for at least some of the observed 
pockets This idea is nominally consistent with Falco’s (1980a) theory concerning pocket formation, 
although the form and motion of the responsible vortex remains in question. 

In figure 23, regions of converging streamlines may be likened to “detachment” lines in a separating 
flow and regions of diverging streamlines to “reattachment” lines. The figure shows that most “detach- 
ment and reattachment” lines lie at a significant angle to the mean (streamwise) velocity vector. Since 
near-wall ±v' motions have a strong association with quasi-streamwise vortical structures (legs), it 

appears that vortex legs are generally oriented at an angle to the mean flow in the x-z plane as well as 
in the x-y plane. 


Pockets with Wall-Pressure Fluctuations 

Figure 24 shows instantaneous streamlines at y+ = 2 together with contours of instantaneous wall 
pressure. As expected, local regions of high pressure coincide with spanwise divergence of the stream- 
lines (wallward motion) and low pressure regions occur along lines of converging streamlines (outward 
motion). (See also Moin and Spalart, 1988.) Thus, pockets may be expected to be evidence of local 
wallward motions above regions of high wall pressure. 
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Bulges in the Outer Interfaces/Large-Scale Motions 

Kovasnay et al. (1970) and others have shown that the instantaneous outer interface between the tur- 
bulent boundary layer and the potential free stream consists of large-scale bulges, separated by deep 
narrow incursions of free-stream flow into the layer. In the simulated boundary layer, the instantaneous 
edge of the rotational flow can be identified by contours of low total vorticity magnitude 

JVcoi + + coij, as shown in figure 25. The picture is essentially the same as given by Spalart (1988), 

and exhibits broad outward bulges and narrow irrotational incursions. The interface bulges have dimen- 
sions of the scale of the boundary layer thickness in both the streamwise (x) and transverse (z) direc- 
tions. Figure 25 resembles visualizations of smoke-filled low Reynolds number boundary layers (e.g. 


The large-scale motion beneath bulges in the simulated boundary layer show a weak rotation in the 
direction of the mean shear, when observed from a frame of reference moving at 0.8 U e . This is in 
qualitative agreement with experimental findings for the large-scale motions (e.g., Blackwelder and 
Kovasznay, 1972; Falco, 1983). 


Backs 


Delta-scale, sloping structures with strong local values of du’/dx have been investigated by a num- 

jdi ? UpS USing multi ' sensor temperature and velocity probes (e.g., Brown and Thomas, 1977; Chen 
and Blackwelder, 1978; Subramanian et al. 1982). Space-time correlations and conditional sampling 
procedures have given from 12 to 30 degrees for the outer-region angle between the wall and these large 
structures. Because these structures are commonly associated with the upstream side of the large-scale 
motions discussed above (Kovasznay et al., 1970), we will refer to them as “backs” (though some publi- 
cations refer to them as “fronts”). 

At this time, we are drawing a distinction between “backs” (which also have transverse vorticity) and 
the previously discussed near- wall shear layers. Backs are tall (delta-scale) discontinuities in u, with 
limited (10-40 viscous lengths) streamwise thickness, and with spacing in the streamwise direction of the 
order of the boundary layer thickness. Near-wall shear layers have dimensions on the scale of the near- 
wall layer, and exist mostly below y + = 80. It is not yet clear if the size and spacing of these shear-layer 
structures are distributed smoothly from small to large, or if there is a bimodal distribution, suggesting 
two different types of structure. Recent statistical analyses of simulated turbulence by Moin et al. (1987) 
ave shown some evidence of a “two-layer” structural makeup of a low Reynolds number channel, 
xperimental ensemble averages have not generally separated the data into two structures possibly 
because the commonly used VITA technique will trigger on either, independent of vertical scale. Also, 
the scale separation between the two structures, if such a separation exists, would be most pronounced at 
high Reynolds numbers, where probe sensors are usually too long to independently detect the small- 
scale structures near the wall. Recent multi-point detection techniques (e.g., Britz and Antonia 1987- 
Fernando et al., 1987) should help clarify this point. 

Following Thomas and Brown (1977), figure 27 shows instantaneous velocity vectors in an x-y 
plane moving at 0.8 U e in the simulated boundary layer. A back structure is clearly evident as a large- 
scale, sloping shear layer with transverse vorticity. The picture is strikingly similar to Thomas and 
Brown s ensemble-averaged velocity field in the vicinity of backs. Figure 28 shows the instantaneous 
transverse vorticity associated with the back in figure 27. 
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Backs with Vortical Structures 


Several researchers have reported on outer-flow vortical structures that are inclined to the wall 
downstream at approximately 45 degrees (e.g., Head and Bandyopadhyay, 1981; Kim and Moin, 1986). 
Hot-wire rake studies have consistently given from 12 to 30 degrees for the slope of the backs of large- 
^J^otions m the outer region (Brown and Thomas, 1977; Rajagopalan and Antonia, 1979; Robinson 
1982). Bandyopadhyay (1980) attempted to resolve this discrepancy by speculating that the shallow- 
angle backs are composed of the heads of arrays of 45 degree hairpin vortices. 

In the current work, an example of the spatial relationship between a back and a hook-like vortical 
structure is shown in figure 26. In this example, the back, which is sloped at about 30 degrees, is located 
upstream of the 45 degree vortical structure. Apparently, the low-speed fluid being left behind by the 
action of the vortex is crowded from behind by the faster upstream flow. The interface between the 
upstream fluid and the low-speed remnants of vortex induction forms the velocity discontinuity which is 
the back. In this scenario, the back is generated without need for the row of hairpin vortices suggested 
by Bandyopadhyay. Final conclusions on the relationships between backs, large-scale motions, and 
vortical structures will require the study of more examples, especially in the simulation at the higher 
Reynolds number, Re 0 = 1410. The generality of the descriptions of the large scale motions in low 
Reynolds number simulations may be limited in view of Murlis et al.’s (1982) finding that the outer 
structure depends on viscosity up to Re 0 = 5000. 


CONCLUDING REMARKS 


The current ongoing investigation of the quasi-coherent structures 
lent boundary layer has produced the following preliminary results: 


in a numerically simulated turbu- 


1 . Examples of all known classes of boundary layer structure have been identified in the simulated 
boundary layer. 


2. High-speed regions in the sublayer are observed to be less elongated than the low-speed regions. 

3. Low-pressure regions tend to be elongated, and apparently serve as a useful spatial detection 
method for instantaneous vortex cores. 


4. Arch-like (head + neck) vortical structures with a broad range of sizes are common in the bound- 
ary layer, while quasi-streamwise vortices (“legs”) with significant streamwise extent are relatively rare. 

, . 5 ‘ High-pressure regions are spheroid and rarely elongated, and typically occur just upstream and 
below the heads of vortical arches, at the interface between high- and low-speed regions. 

6. There is a strong instantaneous spatial association between vortical structures and both elections 
(UV2) and sweeps (UV4). 
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7. Both ejections and sweeps occur along quasi-streamwise “leg” vortices, but ejections are most 
common just upstream and below vortical arches. Sweeps occur most often on the outboard sides of the 

necks” of vortical arches. 

8. Some near-wall shear layers have been observed to roll up into transverse vortices, which are in 
turn associated with significant local contribution to the Reynolds shear stress through ejection (uvo) 

motions. Shear layers not in the process of rolling up contribute only weakly to the Reynolds shear 
strcss< 


9. Backs of outer-region, large-scale motions have been identified in the simulation data. In some 
cases at least, the backs (at a characteristic angle of 18 to 30 degrees) are the interface between low- 
speed fluid trailing a single large-scale vortical structure and the following high-speed fluid. 

Work on the spatio-temporal relationships between the various structures is under way, and results 

the inte™ P t° rted ™? year ‘ ^ m ° ng tbe Questions of primary focus are (1) the nature and importance of 
the interaction of the outer, large-scale motions with turbulence-producing events near the wall; and 

Linen, g ® neratl ° n mech ^ msms and evolution histories of the common vortical structures. As noted, a 
S? P ° P ? a 2° n and , c ° ntnbutl °" statisti « are being compiled to establish the statistical relevance 
Z ; n ! i bS , er !, ed s " features and interactions. Finally, simulated probes and visual markers are 
observa^Foi? t0 m ° re compansons between the simulated boundary layer and experimental 

The current work, along with the 1 987 Center for Turbulence Research Summer Program, marks the 
opening of a new era in experimental fluid dynamics, in which full-simulation results are studied by 
ex^iraenmlists working in “numerical wind tunnels.” The long-term benefits of this new activity 
should include 1) increased experimental emphasis on flows with more complex boundary conditions 
and higher Reynolds numbers, 2) development of advanced instrumentation with specialized measure- 
ment capabilities, 3) increased understanding of probe resolution, conditional averaging, and phase-jitter 
complementary 01 ” 1 expenmenUd/numerical investigations of flows in which the two techniques are 

io l£ 1 L h ° ped , that tbe j current stu dy» as part of a community-wide assessment of the state of boundary 
rf y ^JayCT^uIOT^ 8e, W1 1 aCCelerate progress toward a broad-based understanding of the dynamics 
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FIGURE CAPTIONS 


Figure 1 Computational domain. 

Figure 2- Strategy for structural analysis of numerical simulations. 

Figure 3.- Computational subvolume. 

Figure 4.- Quadrants of the u V plane. 

5.— Contours of fluctuating streamwise velocity in an x-z plane at y - *" = 2. 
Blue to white: u’ + = +0.7 to +2.0. 

Red to yellow: u' + = -0.7 to -2.0. 

Figure 6 - Contours of u’w '2 and u'w ’3 in an x-z plane at y + = 15. 

Green to white: u'w'2/ u| = 2.1 to 14.9. 

Red to white: u'w'3/ = -2. 1 to -14.9. 

Figure 7.- Quadrant 2 shear product in subvolume S. 

Red: uV2/u? = -4.2. 

Figure 8 - Quadrant 4 shear product in subvolume S. 

Blue: u'v'4/ u| = -4.2. 

Figure 9 - Low- speed streaks, ejections, and sweeps in an x-z plane at v + =15 
Yellow: u'+ = -3.0. 

Red: u'v'2/ u| < -4.2 . 

Blue: u'v'4/ u| < -4.2 . 

Figure 10.- Contours of ±v' in an x-z plane at y + = 15. 

Magenta: v’ + = +1.0. 

Green: v’ + = -1.0. 

Figure 11.— Fluctuating vorticity in subvolume S (lowest y"*“ shown is 28). 
Yellow: co^ = ±0. 

Red: 0 )y + = + o. 

Blue: C 0 z + = ±.0. 

Figure 12.- y-z plane cut through low/high-speed streak interface. 

Yellow: u' + = 0 to -b.0. 

Blue: u’+ = 0 to +6.0. 

Green: instantaneous v'w' vectors. 
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Figure 13 - y-z plane cut through low/high-speed streak interface. 

Magenta to white: p ,+ = 0 to -12.0. 

White contour line: p' + = — 4 . 2 . 

Green: instantaneous v'w' vectors. 

Figure 14 - Regions of low pressure in subvolume S. 

White: p'+ = -4.2. 

Figure 15.- Regions of low pressure, low streamwise velocity, and quadrant 2 shear product in 
subvolume S. 

White: p'+ = -4.2. 

Yellow: u’ + = -3.0. 

Red: u'v'2/ = -4.2 . 

Figure 16.- Regions of low pressure, quadrant 2 shear product, and quadrant 4 shear product in 
subvolume S. 

White: p’+ = -4.2. 

Red: u'v'2/ u^ = -4.2 . 

Blue: u'v'4/ u^ = -4.2 . 

Figure 17.- Top view of Figure 16. 

Figure 18 - Near-wall shear layers in an instantaneous x-y plane. 

Blue to green: coj = 0.2 to 1.0. 

Figure 19 - Roll-up of near- wall shear layer. 

Blue to green: co+ = +0.2 to +1.0. 

White: p'+ = -4.2 to -17.0. 

Magenta to yellow: u’v’2/ uf = -2.1 to -12.8 . 

Figure 20.- RoU-up of near- wall shear layer. Instantaneous uv velocity vectors in an x-y plane moving 
X* 

Figure 21.- Regions of high and low pressure in subvolume S. 

White: p’+ = —4.2. 

Green: p ,+ = +4.2. 

Figure 22 - Instantaneous wall-pressure fluctuations. 

Blue to green: p£ = +2.1 to +21.3 . 

Magenta to white: Pw = -2.1 to -21.3 . 


Figure 23.- Instantaneous streamlines in an x-z plane at y+ = 2. 
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Figure 24.- Instantaneous streamlines in an x-z plane at y + = 2, underlaid with contours of wall 
Yellow to red: p ,+ = +3.0 to +25.0. 

Blue to white: p' + = -3.0 to -25.0. 

Figure 25 - Contours of total vorticity magnitude in an x-y plane. 

Blue to magenta: | co + | = 0.02 to 1.0. 

Figure 26.— Back of large-scale motion with regions of low pressure in subvolume S. 

White: p'+ = -4.2. 

Blue to magenta: co z = +0.2 to +3.0. 

Figure 27- Instantaneous uv vectors in an x-y plane moving at 0.8 u e . 

Figure 28 - Instantaneous transverse vorticity in same x-y plane as for Figure 27. 

Blue to magenta: co z = +0.2 to +3.0. 


pressure. 
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Part B - Numerical Simulation Analyses 

Phase 



Strategy 
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Fig. 2. Strategy for structural analysis of turbulence simulations. 



Fig. 3. Computational subvolume. 


2 

EJECTION 


1 

OUTWARD 

INTERACTION 





254 


ORIGINAL PAGE 
COLOR PHOTOGRAPH 





o> o> 


LL LL 



PRECEDING PAGE BLANK NOT FILMED 


29 


„ w ktumHmx mla ** 


4000 







(PUN VIEW) 


ORIGINAL PAGE 
COLOR PHOTOGRAPH 








VECTORS IN CROSS-PLANE 



: low pressure 








PRECEDING PAGE BLANK 


35 


mumhkkt mum* 







PRECEDING PAGE BLANK NOT FILMED 








l U\\w 
1 1 \ 

UH\U\v\ \\\ x i f * 

1 1 % nmtuuiii/^ #, i 
1 1 1 it i imiffifii/** 

1 1 u ittmmmimu^" 
t$titntnfnm$wZTX, 
, 11 / a ttftmmnnnu 
1 1 uuitmmmm\\u.\. 
lit n u utimtmmm,- 
} l % r 
Mil i u 

uut mmmimmt md 
1 1 1 m numimmm^A 
i i inntttMtmtimJ 
* Htiitmmmmn, 3 

XI 1 1 f I f i If fiiftimmm, 
IMmmMimmm* 

Utmmrn* 

{I ij'fttft?**,*' 

; / f ##****, .^ 

I i i t $ i 
ilMtw, 

:}Uii **«,,,, 

f iil n»*,, 4i 

lllln!:::::' 


o 



>% 



\ 


IO <0 

CM CM 


O) O) 


LL LL 


PRECEDING PAGE BLANK NOT FILME 


JUlIEMIMNIMiX §LA»I 







MSA 

National Aeronautcs and 
Space Administration 


Report Documentation Page 


1. Report No. 

NASA TM-102191 

4. Title and Subtitle 


2. Government Accession No. 


3. Recipient's Catalog No. 


A Review of Quasi-Coherent Structures in a Numerically 
Simulated Turbulent Boundary Layer 


5. Report Date 

May 1989 


6. Performing Organization Code 


7. Author(s) 


S. K. Robinson, S. J. Kline (Stanford University, 
Stanford, CA 94305), and P. R. Spalart 


8. Performing Organization Report No. 

A-88294 


9. Performing Organization Name and Address 

Ames Research Center 
Moffett Field, CA 94035 


10. Work Unit No. 

505-60-11 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes ““ ~ — 

Point of Contact: S.K. Robmson^Anws Research^Center, MS 229-1, Moffett Field, CA 94035 

Presented as: Quasi-Coherent Structures in the Turbulent Boundary Layer: Part II. Verification 
and New Information from a Numerically Simulated Flat-Plate Layer at Zoran P. Zaric Memorial 
International Seminar on Near-Wall Turbulence, May 16-20, 1988, Dubrovnik, Yugoslavia 


16. Abstract 


This paper presents preliminary results of a comprehensive study of the structural aspects of 
a numerically simulated number turbulent boundary layer. A direct Navier-Stokes simulation of 
a flat-plate, zero pressure gradient boundary layer at Re 0 = 670 was used. Most of the known 
nonrandom, coherent features of turbulent boundary layers are confirmed in the simulation and 

aS ? eCtS ° f their Spatial character ar e reported. The spatial relationships between 
many of the various structures are described, forming the basis for a more complete kinematical 
of boa^ary layer physics than has been previously known. In particular, the importance 
Qt vortex structures of various forms to the generation of Reynolds shear stress is investigated 


17. Key Words (Suggested by Author(s)) 

Turbulence 
Numerical simulation 
Boundary layer 
Coherent structures 


18. Distribution Statement 

Unclassified-Unlimited 


Subject Category -34 


19. Security Classif. (of this report) 

Unclassified 


— — — 1 ; 

20. Security Classif. (of this page) 

21. No. of pages 

22. Price 

Unclassified 

41 

A03 


NASA FORM 1626 OCT 86 


For sale by the National Technical Information Service Springfield Virginia 22161 




